Influence of Correlated Hybridization on the Conductance of Molecular Transistors 
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We study the spin- 1/2 single-channel Anderson impurity model with correlated (occupancy de- 
pendent) hybridization for molecular transistors using the numerical renormalization-group method. 
Correlated hybridization can induce nonuniversal deviations in the normalized zero-bias conductance 
and, for some parameters, modestly enhance the spin polarization of currents in applied magnetic 
field. Correlated hybridization can also explain a gate-voltage dependence to the Kondo scale similar 
to what has been observed in recent experiments. 
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With the technological advances allowing for the 
construction of single-molecule transistors based upon 
transition-metal complexes [H, [1, there is consider- 
able interest in the out of equilibrium properties of the 
Anderson model. This model is used to describe the cou- 
pling of electrons in localized and strongly interacting 
states of an atom or molecule with electrons in weakly 
interacting extended states of the leads or host metal (see 
the schematic of Fig. 1). Most interest has focused on 
the resonant transport regime of the model associated 
with the local moment or Kondo limit of the Anderson 
model. In this case, correlated (occupancy dependent) 
hybridization has been invoked as a potential explana- 
tion for anomalies in the conductance the unusual 
gate- voltage dependence of the measured Kondo temper- 
ature scale not described by the simplest version of the 
Anderson model, and to suggest a possible conductance 
enhancement via local pairing effects 0|. In parallel, it 
has been demonstrated for transition metal-complexes 
that the correlated hybridization can be very large, com- 
parable to the single-particle hybridization through the 
interatomic potential [6|. 

In this paper, we study the nonequilibrium spin- 
1/2 single-channel Anderson model with correlated hy- 
bridization with the numerical renormalization-group 
(NRG) approach 0, i, i]. We compute the univer- 
sal zero-bias conductance vs. temperature curve of the 
Kondo regime, and show that deviations from univer- 
sality can be induced by correlated hybridization even 
when the magnetic susceptibility appears universal. We 
find, as predicted earlier Q that the spin polarization 
of the currents can be modestly enhanced by correlated 
hybridization in some parameter regimes, though not 
clearly enough to be useful in spintronics applications. 
We also find that the particle-hole symmetry breaking 
associated with correlated hybridization can induce a de- 
pendence of the Kondo scale upon gate voltage similar to 
what is observed in experiment 

To describe the couplings of electrons in the single- 
molecule transistor, we assume the leads are identical 




FIG. 1: Schematic for the spin-1/2 Anderson model with 
correlated (occupancy dependent) hybridization. Here t\ de- 
scribes transfer of an electron into an empty impurity orbital 
from the leads, and t2 describes transfer of an electron into a 
singly occupied impurity orbital from the leads. 



and consider the Hamiltonian for the spin-1/2 Anderson 
impurity model with correlated hybridization, as shown 
schematically in Fig. 1, 

U = ^ek4,<TCk,<T + ^edrfJ-rffT + C/"-dT"di (1) 
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where a conduction electron with momentum k and spin 
a corresponding to an energy eigenvalue ek is created by 
and c^Q creates a conduction electron with spin a 



''k,a' 

at the impurity site 
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The operator d\ creates an electron with spin cr and en- 
ergy Cd at the d level of the impurity, while U is the 
Coulomb interaction, and nd,a ~ d^^da- is the occupancy 
of d orbital with spin a. 
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FIG. 2: (color online) Plots of ksTxiT) / {gusf vs. T /Tk 
for the Anderson model with correlated hybridization. We 
fix Kondo temperature Tk as the temperature for which 
ksTx/igiJ-B)'^ = 0.07. The parameters used are U — 10, D = 
50 and as shown in the legend. 



To solve the Haniiltonian, the numerical renormaliza- 
tion group method provides an excellent approach by 
mapping the Haniiltonian onto a tight-binding Hamilto- 
nian describing a semi-infinite chain with the impurity at 
the end. By using the logarithmic discretization and di- 
agonalizing the chain Hamiltonian iteratively while keep- 
ing the lowest eigenstatcs, one can calculate the thcrmo- 
dynamical and dynamical properties with good precision. 
As a result of the logarithmic discretization, the hopping 
along the chain decreases exponentially; the hopping am- 
plitude between site n and site n -I- 1 is proportional to 
A^"/^, where A is the discretization parameter. 

In the quest for new fixed points of the model we inves- 
tigated the NRG level flow as well as the temperature- 
dependent static magnetic susceptibility. In our calcula- 
tion, we set the discretization parameter A = 1.65, and 
keep the lowest 1500 states in each iteration. The band- 
width D is 50, and the Coulomb repulsion U is set to 
10. The energy scale is fixed by setting the Fermi- level 
density of states per spin pp = I/tt. 

Figure 2 shows the value of ksTxiT) / (gfJ-B)^ as a func- 
tion of T/Tk- In the presence of correlated hybridiza- 
tion, the drop of Tx from the local-moment regime or 
mixed- valence regime to the strong-coupling regime looks 
exactly the same as the temperature decreases. By an- 
alyzing the level flow, we find that there are no novel 
fixed points induced by correlated hybridization; all the 
effects are rolled into a suitable renormalization of the 
Kondo scale. The correlated hybridyzation may, how- 
ever, change the steepness of the drop of Tx if the system 
goes from the free-orbital regime directly to the frozen- 
impurity regime, such as the cases of ti = t2 = 1, = — 1 



FIG. 3: (color online) Normalized conductance G{T, V = 
0)/Go, where Go is the zero-temperature conductance, as a 
function of T/Tk- 



and ti = 0.7, t2 — 3, ej, — ~3, where there is no Kondo 
universality. 

The conductance for the single-molecule transistor 
with symmetric leads at zero bias can be expressed as 
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where pf is the Fermi- level density of states, /(w) is the 
Fermi distribution function, and Pa{^) is the transimis- 
sion density spectral function. We carry out the calcu- 
lation of transimission density Po-(^) in Lehmann repre- 
sentation. 
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where the operator A„ of the Green's function 
{{AMl)){^ ~ iS) is given by A„^h\Q){cT\- t2\-a){2\ 
with |0), |(t), and |2) being the empty, singly, and doubly 
occupied impurity states, respectively. Z is the partition 
function, and |n),|m) are the eigenstates of the Hamil- 
tonian in Eq. ([1]). Therefore, the conductance can be 
obtained by 



(6) 



Since Pai^^) only contributes in the Fermi window, the 
high energy excitation do not contribute, and it is good 
to take into account only the states of the last iteration 
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of the NRG calculation. The influence of correlated hy- 
bridization on the conductance is investigated by study- 
ing the Cd dependence of the conductance. In the molec- 
ular transistor, ed can be tuned with gate voltage. As 
a consequence of correlated hybridization, the symmetry 
of the <-> 1 and 1^2 valence fluctuations in the con- 
ductance is broken, with the resonance widths propor- 
tional to accompanied by the shifts of the resonance 
positions. As the temperature decreases from high tem- 
perature to the Kondo temperature, we flnd the enhance- 
ment of the conductance, the broadening of the conduc- 
tance peaks, and the shifts in their positions towards each 
other, which are characteristics of the Kondo effect. In 
the limit of T ^ 0, the conductance per spin calculated 
by the NRG method is found to satisfy the unitary limit, 

2 o 

Ga = j;^sm Sa, and the peak of total conductance is 
always found Gmax = 2e^//i to within A-dependent sys- 
tematic corrections for different magnitudes of correlated 
hybridization. In the other words, the phase shift analy- 
sis using Friedel sum rule holds with correlated hy- 
bridization; the Fermi energy phase shift for an electron 
of spin (T scattering off the impurity site can be obtained 
by Scr = ■n{nd,a)- 

It has been believed that the conductance normalized 
to its zero-temperature value is universal in the Kondo 
regime [13, III as well. Figure 3 shows G/Gq, where 
Go is the zero-temperature conductance, as a function of 
T/Tk- We flnd that for parameters within the range of 
Kondo regime, such as ii = 0.3, ^2 = 1, — 5 < < —1, 
the curves fall on top of those without correlated hy- 
bridization, showing the Kondo universality. However, 
we notice that in the case of ti = 0.7 and <2 = 3, which 
has a magnitude of \a) — > |2) transition comparable to U, 
the conductance curves are obviously off the universality, 
while there is only slight variation in the Tx curves, as 
shown in Fig. 2. We also note that a pcrturbativc treat- 
ment of correlated hybridization in Anderson Hamilto- 
nian ^ was used to explain the anomalous conductance 
plateau around 0.7(2e^//i) observed in experiments on 
quantum point contacts. We do not see the 0.7 anomaly 
in the conductance calculated by this nonperturbative 
approach. In addition, there is no indication of the en- 
hancement of the conductance in our calculation induced 
by local pairing Q associated with correlated hybridiza- 
tion. 

Correlated hybridization has been proposed using 
mean-field treatment as a potential means to spin polar- 
ization of currents through transition-metal based molec- 
ular transistors in modest magcnctic fields However, 
the effect on spin polarization in the low-temperature 
limit is still not clear. Here, we study the ratio of spin 
conductance at small magnetic field to total conductance 
at zero field in the zero-temperature limit, ^Q^j^^Q-j , as a 
function of single-particle energy. The magnetic field is 
set to H = Q-lksTK / {gus) for each set of parameters. 
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FIG. 4: (color online) The ratio of spin conductance at finite 
field, H = GAksTK / {gyiB), to total conductance at zero field 
in the zero-temperature hmit as a function of single-particle 
energy. 



as shown in Fig. 4. Without correlated hybridization, 
i.e. <i = t2 = 1, the ratio is zero at = —5, where 
there is particle-hole symmetry, and the zero-field phase 



shift (5o = 7r/2. The ratio 



increases to about 0.06 



G(H=Q) 

and saturates in the empty and doubly occupied regime. 
When ti is reduced, the position of (5o = 7r/2 is shifted 
towards e^; = because the larger hybridization between 
singly and doubly occupied impurity states renormalizes 
their energy levels down relative to the level of empty 
impurity state. The spin polarization of currents can 
be enhanced depending upon the value of e^. When 
ti = 0.7 and t2 = 3, there is noticeable enhancement in 
the spin polarization in the regime between = —3 and 
Ed = —8. Interestingly, the saturated value of the ratio of 
spin conductance to zero-field total conductance remains 
the same for different magnitudes of hybridization, and 
it seems to have universality. 

One of the most interesting results with correlated hy- 
bridization is how it affects the dependence of the 
Kondo scale. Generally without correlated hybridiza- 
tion, as shown in Fig. 5(a), the Kondo temperature 
decreases when moves away from the charge degen- 
eracy points, Erf ~ and ed = —10, towards the particle- 
hole symmetry point, = —5. In the limit of F <C J7, 
the Kondo temperature has been well described by scal- 
ing theory [H, Tk ^ l/2Vr!7exp[7red(ed + C/)/(rC/)], 
in the Kondo regime. In the presence of correlated hy- 
bridization, the position of the minimum of Kondo tem- 
perature shifts away from Cd ~ —5. The shift is espe- 
cially significant in the case of ti = 0.7 and t2 = 3, 
where there is less dependence of Kondo tempera- 
ture. The particle-hole symmetry breaking associated 
with correlated hybridization may explain the anoma- 
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ti = 0.3, ^2 = 1 when e^/F moves away from 0. The de- 
crease of Tk becomes less dramatic when e^/r goes into 
the Kondo regime (e^/r ^ —1), and the dependence be- 
comes much weaker when e^/r < —2. The curve seems 
to capture the feature of the weak gate-voltage depen- 
dence of the Kondo temperature in the single-molecule 
transistor found by Liang et al. [H We note that the or- 
bital degeneracy which plays a role in Ref. [H is absent 
in our model, so our calculation of Kondo temperature 
would have more dependency on e^/r. 

The present study shows that the correlated hybridiza- 
tion using spin-1/2 single-channel Anderson model can 
induce deviations in the normalized zero-bias conduc- 
tance from universality, and it can modestly enhance the 
spin polarization of currents in small magnetic field. The 
particle-hole symmetry breaking associated with corre- 
lated hybridization provides an explanation for the weak 
gate-voltage dependence of the Kondo scale similar to 
what has been observed in transistion-metal based molec- 
ular transistors. Further investigations of correlated hy- 
bridization in the degenerate Anderson model are ex- 
pected to provide detailed insight into electron transport 
through a more realistic singlc-nioleculc transistor with 
orbital degeneracy. 
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FIG. 5; (color online) (a) Kondo temperature as a function 
of single-particle energy, (b) Tk (normalized) as a function 
of ed normalized by the average hybridization F. The legend 
is the same as that in (a). 

lous results in a recent Letter Q that the Kondo tem- 
perature is less dependent on gate voltage and even in- 
creases when moving away from the supposed charge 
degeneracy point, = 0, in single- molecule transis- 
tors based upon transition-metal complexes. In fact, 
the true charge degeneracy point is shifted away due to 
the renormalization of the energy levels of the impurity 
states in the presence of correlated hybridization. Fig- 
ure 5(b) shows rft:/rA'(e(i/F=— 1) as a function of e^/F, 
where F = TTp[ti{l ~ {rid^a)) +t2{nd,iy)^ is the average 
hybridization which has occupancy dependency. How 
the Kondo scale changes with e^/F depends upon the 
magnitude of correlated hybridization. We see that in 
the mixed- valence regime (— 1 < Cd/F < 0) the Kondo 
temperature decreases by about six folds in the case of 
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